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Abstract 

We consider the effective lagrangian describing the J\f = 1 supersymmetric 
QCD with SU(N C ) gauge group involving Nf = N c — 1 pairs of chiral mat- 
ter multiplets in fundamental and antifundamental color representations. For 
this theory in the framework of the effective lagrangian approach, we solve 
BPS equations for the domain walls interpolating between different vacua. 
The equations always have a unique solution for the walls interpolating be- 
tween the chirally asymmetric and a chirally symmetric vacua. For the walls 
interpolating between different chirally asymmetric vacua, the equations ad- 
mit two different solutions which exist in a limited range of the mass of the 
matter fields m < to* = kAsqcd where the parameter k depends on N c . At 
m = to*, two branches join together and, at to > to*, no BPS - saturated 
complex domain walls exist. 



1 Introduction 



Supersymmetric QCD is the theory involving a gauge vector supermultiplet V and 
an even number of chiral matter supermultiplets Si, S[ (i = 1, . . . , Nf is the flavour 
index) such that Si belonging to the fundamental representation of the gauge group 
SU(N C ) involve left quarks and their scalar superpartners and S"- belonging to the 
antifundamental representation involve left antiquarks and their scalar superpart- 
ners. The lagrangian of the model reads 

C = ^Tr J d 2 6 W 2 + H.c.^ + EQ/ d20d ^ Sie v St 

+ \J #m S[e~ v S[) - / d 2 9 S' t S, + H.c.) , (1.1) 

rriij is the complex mass matrix; color and Lorentz indices are suppressed. 

The dynamics of this model attracted the attention of theorists since the be- 
ginning of the eighties. In some respects, it is similar to the dynamics of the stan- 
dard (non-supersymmetric) QCD. If the number of flavours Nf is not too large, 
it involves confinement, for instance. The model (|1.1|) has, however, many specific 
features which are due to supersymmetry. Actually, this theory is in some sense 
much simpler than the standard QCD — the presence of the extra symmetry allows 
one to obtain many exact results. The results are especially strong for the extended 
M = 2 version of the theory but a lot of exact theorems can be derived also for 
M = 1 theory @. 

The details of the dynamics depend on Nf and on the matter mass matrix m^-. 
A common feature of all models (|1.1|) is their nontrivial vacuum structure associated 
with the spontaneous breaking of a discrete chiral symmetry. Like in the standard 
QCD, the axial C/a(1) symmetry corresponding to the chiral rotation of the gluino 
field and present in the tree-level lagrangian ( |1 . 1| ) is broken by anomaly down to 
Z2N C - This discrete chiral symmetry can be further broken spontaneously down to 
Z 2 so that the chiral condensate < Tr {\ a \ a } > (A Q is a two - component Weyl 
spinor describing the gluino field) is formed. There are N c different vacua with 
different phases of the condensate 

< Tr A 2 > = £e 2 " fc/A S k = 0, . . . , N c - 1 (1.2) 

It was noted recently || that on top of N c chirally asymmetric vacua (|1.2| ), also a 
chirally symmetric vacuum with zero value of the condensate exists. 

The presence of different degenerate physical vacua in the theory implies the 
existence of domain walls — static field configurations depending only on one spatial 
coordinate (z) which interpolate between one of the vacua at z = — oo and another 
one at z = oo and minimizing the energy functional. As was shown in in many 
cases the energy density of these walls can be found exactly due to the fact that the 
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walls present the BPS-saturated states. []. 

One can show (see Refs.|5], || for detailed derivation and discussion) that the 
energy density of a BPS-saturated wall in SQCD satisfies a relation 
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where the subscript ±00 marks the values of the gluino condensate at spatial infini- 



ties. The RHS of Eq.(L3) presents an absolute lower bound for the energy of any 
field configuration interpolating between different vacua. 

The relation (|1.3| ) is valid assuming that the wall is BPS-saturated. However, 
whether such a BPS-saturated domain wall exists or not is a non-trivial dynamic 
question which can be answered only in a specific study of a particular theory in 
interest . 

In Refs.l 



TO, ITTJ this question was studied in the simplest nontrivial case iV c 



2, Nf = 1 in the framework of the effective low energy lagrangian due to Taylor, 



Veneziano, and Yankielowicz []12j] . In that case, it is written in terms of the composite 
colorless chiral superfields 
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and presents a Wess-Zumino model with the superpotential 



W = -$ 3 
3 



The results are the following: 
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1. For any value of the mass of the matter fields m, there are domain walls 
interpolating between a chirally asymmetric and the chirally symmetric vacua. 
They are BPS - saturated. 

2. There are also complex BPS solutions interpolating between different chirally 
asymmetric vacua. But they exist only if the mass is small enough m < m* = 
4.67059 . . .Asqcd- When m > m*, BPS walls are absent. 

3. In a narrow range of masses < m < m„ m A.83Asqcd, complex domain 
walls still exist, but they are not BPS saturated anymore. At m > m**, there 
are no such walls whatsoever. 

All these results were obtained by solving numerically the first order BPS equations 



e l5 dW/d<f), d zX 



;i.6) 



1 Such BPS-saturated walls were known earlier in 2-dimensional supersymmetric theories (they 
are just solitons there) || (7) and were considered also in 4-dimensional (non-gauge) theories in 
stringy context M. 
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associated with the TVY lagrangian Q (see Refs.[|, [| for details) and/or the equa- 
tions of motion for the profile of the wall with the proper boundary conditions. 
Technically, the easiest problem was to study the point 1 of the list above ||: the 
phase of the complex scalar fields (f) and \ (the lowest components of the superfields 
$ and X) does not change along such a wall and can be set to zero which simplifies 
the equations considerably. Point 2 is a bit more difficult [flO|j : one has to solve here 
essentially complex BPS equations to obtain essentially complex solutions. Point 
3 is still more difficult: we have to solve the equations of motion which are of the 
second order [U]] . 

In this paper, we address the same issues for theories with higher unitary gauge 
groups. In particular, we study the points 1 and 2 and show that the results are basi- 
cally the same: the real walls are always BPS saturated, and they exist for all masses 
while the complex BPS equations have non-trivial solutions only in the limited range 
of mass: m < m* = .28604. . .Asqcd for SU(3), m < m* = .07539 . . -Asqcd for 
SU(4), etc. When m > m*, the solutions are absent. A study of the corresponding 
equations of motion which would allow us to find the form and the properties of 
non-BPS walls is now in progress. 

The paper is organized as follows. In the next section we write down the effective 
lagrangian for the model in interest, discuss its vacuum structure and some general- 
ities concerning the wall solutions. In Sect. 3, we find the BPS solutions for the real 
walls. Complex walls in the Higgs phase (in the limit m — > 0) are discussed in Sect. 
4. In Sect. 5, we analyze the case of small finite masses and find the solution as a 
series in a small Born-Oppenheimer parameter. Numerical results for the complex 
walls for general values of mass are presented in Sect. 6. In Sect. 7, we discuss the 
properties of a peculiar "lower" BPS solution at small masses. Summary and some 
general discussion is the subject of the last section. 



2 The model 

In this paper, we restrict ourselves with studying the model ( |1 . 1| ) with Nf = N c — 1. 
We also choose the simplest form for the mass matrix rrtij = m5ij. When m is small 
(much smaller than A S q C d), the expectation values of the squark fields (sj), (s'j) 
in chirally asymmetric vacua are large. It turns out that, for different flavours, the 
vacuum expectation values have different colour orientations. Up to overall flavour 
and colour rotations, one may choose 

\ 
1 

o / 

,1,0) (2.1) 
"direction" of the wall. 
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2 In the SU(2) case, we should set S = or S = it depending on the 
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with v oc m 1 / 27Vc ; columns and rows display the colour structure. A set of Higgs 



averages (|2.1|) break the SU(N C ) colour symmetry completely: all the gauge bosons 
acquire a large mass. When m is small, v is large, the effective coupling g{y) is 
small and the vacuum state has a trivial perturbative nature. One can say that the 
theory is in a weak-coupling Higgs phase. 

The effective TVY lagrangian is written now in terms of the composite colorless 
fields 

3 
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It is again a Wess-Zumino model with the superpotential 
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The superpotential is rigidly fixed from the requirement that the conformal and 
the chiral anomaly of the theory ( |I . 1 1 ) under consideration are reproduced correctly. 
For all vacuum states, the matrix (s^Sj) of the vacuum expectation values is propor- 
tional to unity: (s^Sj) = CSij. The same is true for the domain walls interpolating 
between different vacua. So, we can safely impose the constraint M.^ = X 2 Sij 
and study the Wess-Zumino model with just two superfields <3? and X and the 
superpotential 
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ln^X 2 ^- 1 )} - ll --(iV-l)X 
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(JV = N c and we have set Asqcd = 1 )■ The corresponding potential for the lowest 
components 0, x °f f ne superfields $, X is 
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The potential ( |2.5| ) has N + 1 degenerate minima. One of them is chirally 
symmetric = x = (one should take care, of course, that the ratio 3 /x would 
also tend to zero in the limit <fi — ► 0, x ~ * 0)- There are also A/" chirally asymmetric 



vacua 



with 
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The values of the superpotential (|2.4j) at the minima are W sym = and 




(2.6) 
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There are also the minima with the inverse sign of x (and the appropriately chosen 
phase of 0), but they are physically the same as the minima ( |2.6|) : the vacuum 
expectation values of the gauge invariant operators < Tr A 2 > = (327r 2 /3) < 3 > 
and < 2s-Sj > = < x 2 > $ij are the same. The vacuum values of the gluino 
condensate are (Tr A 2 ) sym = and 



(Tr A 2 ) fc 



4m 



N-l\ 1/N 



in 



D 2mk/N 



(2.9) 



To study the domain wall configurations, we should add to the potential ( |2.5p 
the kinetic term which we choose in the simplest possible form 



C 



kin 



+ m 



(2.10) 



(it follows from the term 1/4/ d 4 6<&& + 1/4/ d A 6XX in the superlagrangian) and 
solve the equations of motion with appropriate boundary conditions. Following 
Refs.p, |10[, we will look in this paper only for the solutions of the more simple BPS 
equations (|1.6|) which are of the first order. 

The model involves two types of walls. Some of them interpolate between the 
chirally symmetric vacuum and a chirally asymmetric one. There are N such walls 
which transform into each other inder trivial Zn phase rotation. It is convenient to 
choose the chirally asymmetric vacuum with k = in which case the wall solution 
is purely real. 



Other walls interpolate between chirally asymmetric vacua (2J:) with different 
k. We will restrict ourselves with studying the walls interpolating between adjacent 
asymmetric vacua, say, between the vacuum with k = and the vacuum with k — 1. 
For N = 2, 3, it is the whole story but, starting from N = 4, also the walls connecting 
the vacua with k = and k = 2 etc. may appear. All such solutions are essentially 
complex. 

The value of 5 to be chosen in Eq. ( |1.6[ ) depends on the wall solution we are going 
to find. To fix it, note that the equations (|1.6|) admit an integral of motion: 



Im[W(0,x)e 
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Indeed, we have 



e l5 d z W = e 15 [ -K-d z( p + -K-d zX 
V d(f) dx 
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The real wall connects the vacua with W = and W = — A^/2(4m 7V_1 /3) 1 / Ar . These 
boundary conditions are consistent with Eq.( |2.1lQ only if e lS = ±1 (the sign depends 
on whether the walls goes from the symmetric vacuum to the asymmetric one when 
z goes from — oo to +oo or the other way round). 
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A complex wall connects the vacua where the superpotential ( |2.4| ) acquires non- 
zero values with different complex phases. For the wall connecting the vacua with 
k = and k = 1, the condition ( |2.11| ) is consistent with the choice 5 = n/N ± n/2 
depending on the direction of the wall. 

Bearing in mind Eqs. ( |1.3j , |2.9| ), the energy densities of the BPS walls are 



(Am N - l \ 1/N , 
N[——\ (2.12) 



for the real walls and 



7T (Am N ~ l \ 1 IT 

e c = 2iVsin — = 2e r sin— (2.13) 

N \ 3 J N v ' 

for the complex walls. 

Before coming to grips with finding the numerical solutions of the BPS equations 
for an arbitrary value of mass, let us discuss what happens in the limiting cases 
m — > oo and m —>■ where the situation is considerably simplified. 

Consider first the case of large masses. In this case, one can integrate the heavy 
matter fields out and write the effective lagrangian for the composite chiral superfield 
$. Technically, one should use the Born-Oppenheimer procedure and to freeze down 
the matter field x so that the large second term in the potential (|2.5| ) disappear. 
Proceeding in supersymmetric way, we get X 2 = 4$ 3 /3m. Substituting it in the 
first term, we obtain the Veneziano-Yankielowicz effective lagrangian [TBI] which is 



the Wess-Zumino model for the single chiral superfield $ with the superpotential 

2N 



W = — $ 3 
3 

where $ is measured now in the units os 



ln<r- 1 



(2.14) 



A SYM - \ ~J A SQCD 

The corresponding potential is 

U{<j>) = \dW/d(j)\ 2 = 4iV 2 |0 4 ||ln0 3 | 2 (2.15) 

This expression is not yet well defined: the logarithm has many sheets, and one 
should specify first what particular sheet should be taken. An accurate analysis 
taking into account the fact that the topological charge v oc / r Ti{G 2 G 2 }d A x in 
the original theory ( |1.1|) is quantized to be integer reveals that the true potential 
is glued out of iV such sheets P, ||. The gluing occurs when the phase of the 

Quite an analogous situation holds in the Schwinger model: the true bosonized lagrangian (in 
the Schwinger model it is just equivalent to the original theory ) is glued out of several branches 
when taking into account the effects due to quantization of topological charge H] . 
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expression under logarithm reaches the values ±7r /N HQ: 
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(remind that only the field <j) 3 oc Tr A 2 , not <fi itself has a direct physical meaning) 



Bearing the prescription (|2.16|) in mind, the potential (|2.15 ) has A^ + 1 minima as 
it should: sym = and fc = e 2i7rk / 3N . 

As was shown in Ref. || , the BPS equations involve only the real solutions with 
the chirally symmetric vacuum at one of the boundaries in this case. Non-trivial 
complex domain walls connecting different chirally asymmetric vacua are absent. 

Consider now the case of small masses. In this case, chirally asymmetric vacua 
are characterized by large expectation values of the matter scalar field (%) ~ oc 
m -i/2N _ Again, the theory involves two different energy scales, and one can tenta- 
tively integrate out heavy fields and to write down the Wilsonean effective lagrangian 
describing only light degrees of freedom. Proceeding in the Born-Oppenheimer 
spirit, we should freeze now the heavy field <j) m the potential fl2.5|) so that the large 
first term in the potential acquire its minimum (zero) value. In contrast to the large 
mass situation, this can now be achieved in two ways: either by setting <p = or 
by setting (f) 3 x 2<yN ~ 1 ^ = 1. In the first case, we will obtain the effective lagrangian 
describing the dynamics of the chirally symmetric phase which is just the lagrangian 
of free light chiral field X involving x an d its superpartner. 

The second choice results in the lagrangian describing the dynamics of the chi- 
rally asymmetric phases. It is the lagrangian of the Wess-Zumino model with a 
single chiral superfield X and a non-trivial superpotential 



W = - 



(2.17) 



This lagrangian is well known and was obtained earlier from instanton and/or from 
holomorphy considerations f2j. The corresponding potential U = \ dW/dx\ 2 has 
different non-trivial minima at Xk = p*e mk ^ N . When m < 1, a large expectation 
value (x) results in breaking down the gauge symmetry of the original theory by the 
Higgs mechanism: the theory is in the Higgs phase. 

The BPS equations corresponding to the superpotential (|2.17|) admit non-trivial 
complex domain wall solutions connecting different asymmetric vacua. When N = 2, 
the solution can be found analytically 0. For A^ > 3, there are only numerical 
solutions which will be presented and discussed in Sect. 4 of the paper. 
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In this approach, we are not able, however, to study the real domain walls 
interpolating between the chirally symmetric and a chirally asymmetric vacua. Such 
a wall corresponds to going through a high energy barrier separating two kind of 
vacua. It is a remarkable consequence of supersymmetry and of the related BPS 
condition ( |1.3| ) that the energy of such a wall is still not large. For N = 2, such 
walls were analyzed in Ref. 0. In the next section, we do the same for arbitrary iV 
and reach the same conclusions. 

Thereby, the physical situation in the limits m — > and m — > oo is somewhat 
different. In both cases, we have N + 1 different vacua. However, an analog of 
nontrivial complex walls connecting different asymmetric vacua which are present 
at small masses, is absent when the mass is large. In the latter case, only real walls 
are present. It is therefore very interesting to understand what happens in between, 
at intermediate values of masses, and how the transition from one regime to another 
occurs. That was the main motivation for our study. 

All the calculations were performed in the framework of the effective theory with 
the potential ( p.5|) . The status of this effective theory is somewhat more uncertain 
than that of ( |2.17| ) — for general value of mass, the TVY effective lagrangian is not 
Wilsonean; light and heavy degrees of freedom are not nicely separated. But it pos- 
sesses all the relevant symmetries of the original theory and satisfies the anomalous 
Ward identities for correlators at zero momenta. We think that the use of the TVY 
lagrangian is justified as far as the vacuum structure of the theory is concerned. 



3 Real Walls 

The BPS equations ( |1.6|) with the superpotential Q2.4|) have the form 



e^^ln^x 2 ^} 
si /AT i V -mx 



X' = e ld -(N-l 



3X 



(3.1] 



(O' = d z O). To find the wall interpolating between = x = Oatz = — oo and 
= i?*, x — P* a t z = oo, we have to choose 5 = n (or 5 = for the wall going 
in the opposite direction). With this choice and the boundary conditions given, the 
solutions <j>(z) and x( z ) are g° m g to be real so that we have a simple system of just 
two first-order differential equations. 

For all N, the dynamics of this system is quite similar to that in the case N = 2 
studied in Ref. ||. The solution exists for all masses. For large m, the heavy matter 
field can be integrated out, and we arrive at the BPS equation for the supersym- 
metric gluodynamics 

4>' = -2N<p 2 \nU 3 /Rl) (3.2) 



The solution of this equation with the boundary conditions (f>(— oo) = 0, 0(oo) = 
can be expressed into integral logarithms. 
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As was the case for N = 2, the solution of the system ( |3.1| ) can be found in 
an analytic form also for small masses. The wall trajectory consists then in two 
distinct regions. On the first stage, only the light field x is changed, x( z ) = P*e mz , 
while the field stays frozen at zero. Then, at z — 0, the trajectory abruptly turns. 
x(z) stays frozen at its vacuum value p* and the equation for <p(z) is reduced to 
0' = — 20 2 ln^ 3 /? 2 /^ -1 ^} which is the same as Eq. ( |3.2| ). The second stretch is 
much thinner than the first one, its width being of order i?" 1 oc m~ ( - N ~ 1 ^ 3N . It is on 
this second stage when the high and narrow potential barrier between the regions 
~ and ~ ^~ 2 ( Ar ~ 1 )/ 3 is penetrated. The width of the first stretch is ~ 1/m 
and it carries the fraction (N — 1)/N of the total energy (|2.12j ). Correspodningly, 



the second thin stretch carries the fraction 1/N of the total energy. 

In the intermediate range of masses, the solution has to be found numerically. 
Parametric plots in the ((f), x) plane for N = 3 and different values of m are drawn 
in Fig. IJ 

4 Domain Walls in Higgs Phase 

We have got to solve the BPS equations 

X' = e-^V 2 ^, (4.1) 
ox 

where W is the superpotential ( |2.17| ), with the boundary conditions 

x(-oo) = p*; x(oo) = p^' N ; (4.2) 

It is convenient to introduce polar variables x = pe m . The equations (|4. 1| ) acquire 
the form 

P = (iV-l)|mpsin 7 - — A-jsin^JV-l)]! 
i = 2(iV-l)|mcos 7 + ^cos[ 7 (iV-l)]| (4.3) 

where 7 = 2a — tt/N changes from 7 = —tt/N at z = —00 to 7 = ir/N at z — 00. 
For N = 2, these equations were solved in Ref.||. The solution is analytic: 

p{z) = p* 

tan j(z) = sinh[4m(z — zq)] (4.4) 



or in the complex form: 



1 + ie 4m(z-z ) 
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Figure 1: Real walls (N = 3). The parametric plots for m — 1 (solid line), m = .1 
(dashed line), and m — 10 (dotted line). and x are measured in units of i?* and 
p*, respectively. 
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Figure 2: BPS walls in Higgs phase for N = 3 (dashed line), N = 5 (dotted line) 
and N — 10 (solid line). 
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(z is the position of the wall center). For N > 3, the solutions can be found 
numerically. The profiles for the ratio r(z) = p(z)/p* in the interval z = < z < oo 
[it is a half of the wall, another half being restored by symmetry considerations: 
p(—z) = p(z)) with different N = 3, 5, 10 are presented in Fig. ||. 

We see that the dependence p(z) is not flat anymore but displays a bump in the 
middle. To understand it, remind that the system ( [4.3| ) has the integral of motion 
( p. 11| ). In our case, it amounts to 

m(N-l) 2 2 r /lir N1 N fAm N - 1 \ 1/N n , . 

2 P cos 7 ^ ^rrycos[ 7 (iV-l)] = _ (^^— J cos- (4.6) 

as follows from the boundary conditions ( f4.2|) and the relations ( |2.8|) . In the middle 
of the wall, 7 = 0, and the condition (|4.6| ) implies 



(AT-1)^__L_ = iVcos^ (4.7) 

( x = r(0)). It is not difficult to observe that the real root of the algebraic equation 
( |4.7D is slightly greater than 1 for N > 3. When N is large, x — 1 tends to zero 
oc 1/N. 



5 Born — Oppenheimer Expansion at Small Masses 

When the gauge degrees of freedom associated with the superfield $ do not 

decouple completely and should be taken into account. The full system (|3.1| ) of the 
BPS equations for the complex domain walls has the form 




a' 

R' 

& 



(N-l) 
(N-l) 
-2R 2 
2R 



mpsm(2a - ix/N) - ^- sin(3/3 - ix/N) 



3p 



mcos(2« - ir/N) - cos(3/? - vr/iV) 



sin(3/3 - tt/N) \n(R 3 p 2 ^ N ^) + cos(3/3 - vr/iV)[3/3 + 2a(N - 1) 
cos(3/3 - ir/N) \n(R 3 p 2 ( N -V) + sin(3,5 - tt/N)[3(3 + 2a(N - 1)] 



(5.1; 



where, as in the previous section, we have chosen 5 = tt/N — it/ 2 and introduced 



the polar variables \ 
boundary conditions 



pe l 



Re l/3 . One should solve the system ( |Q| ) with the 



-oo 



p(-oo) 
(3(-oo) = 0; 



p(oo) 



a oo 



■ p*; R(-oo) 
tt/N; (3(oo) 



R(oo) ■ 
2(N 



1)tt 



3N 



(5.2) 



When N = 2, the system (|5.1|) is reduced to that studied in Refs.[0, p4| . 

Generally, we should solve the full system ( |5.1|) numerically. The problem is 
simplified when noting that the wall solution should be symmetric with respect to 
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its center. Let us seek for the solution centered at z = so that 



p{z) = p(-z), R(z) = R(-z), 
a{z) = n/N - a{-z), (3(z) = -2{N - l)n/(3N) - /3(-z) (5.3) 



Indeed, one can be easily convinced that the Ansatz (|5.3|) goes through the equations 
( |5.1|) . It is convenient to solve the equations ( |5.1| ) numerically on the half-interval 



from z = to z = oo. The symmetry (|5.3|) dictates a(0) = n/(2N), (3(0) = 
— (N — l)n/(3N). The values p(0) and -R(O) are related by the condition ( |2.11| ) 
which reads here 

^^{ H R*(0)p^\0)}-l}+mp*(0)(N-l) = N^^y N cos^(5A) 

Thus, only one parameter at z = [say, R(0)] is left free. We should fit it so that the 
solution would approach the values specified in Eq.( |5.2| ) at z — > oo. The numerical 
solution for this problem will be presented and discussed in the next section. Here 
we will concentrate on the case when m is small where some analytic results can be 
obtained. 

In the limit m — > 0, we can just freeze the heavy variables: 

(3 = -2a(N - l)/3; R = p- 2 ^-^ 3 (5.5) 



in which case the first two equations in Eq. ( |5.1| ) reproduce the system ( fO|) studied 



in the previous section. This is the leading order of the Born-Oppenheimer expan- 
sion. To proceed further, we should allow for the fast variables R, j3 to deviate 
from their zero order values ( ]5.5[) . Integrating these deviations out, one can obtain 
the corrections to the leading order Born-Oppenheimer hamiltonian for the slow 
variables p, a. The spectrum of the low-energy effective hamiltonian and all other 
quantities of interest can be expanded in a series over a small parameter. In our 
case, the relevant expansion parameter is ~ m/R* oc m ( 2N + 1 )/ 3N . 

Such an expansion can be carried out also on the level of dynamic equations. 
One should present 

P = po+Pi+P2+- ■ ■ , a = ao+ai+«2+- • • , R = Ro+Ri+R2+- • • , (3 = Po+Pi+P2+- ■ 



, where p and a are the solutions of the system (|4.3|) and R , /3 are found from 
Eq.([5.5|), and linearize the equations in each subsequent order. Let us do it for 
the BPS system ( |5.1|) . We will be particularly interested in the quantity R(0) to 
compare it with the numerical results of the next section. 

Let us look first at the two last equations in Eq.(|5.1|). In the first non-trivial 
order, we have 



R' = -2Rl {sin(3/? - n/N) 
(3' = 2R {-cos(3(3 -tt/N) 



3Ri , 2(jV-l)pi 

Ro Po 

3-Ri , 2(AT-l)pi 

^0 PO 



+ cos(3/5 - ir/N)[3/3x + 2 ai (iV - 1)]} 
+ sin(3/5 - tt/N)[3(3i + 2 ai (N - 1)]} 



(5.6) 
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From this, one can readily express R\ and (3\ to be substituted in the linearized 
version of the first pair of the equations in Eq. (|5.1|) . As a result, we obtain the 
linear system for p\ and 71 = 2«x — n/N : 



Pi 



(N-l) 



A(N — 1) _2(2JV+l)/3 



9 



Po 



Po + Pi 



4(2iV - 1) 

msin 7oH — sin[7 (A^- 1)] 

6 Po 



7i 



mp COS70 



A(N-l) 

(27V-1) - cos[ 7o (iV-l)] 



7i 



2(N-1] 



3p ' 
2{N - 1 



-7i 



m sin 70 + 



9 

4(JV 



■Po" 2(2 " +1)/3 7o 



8N 



o (2AT+1) 

6 Po 



cos[7o(A^ - l)]pi 



3pg 



2N 



■sin[7o(JV-l)] 



The proper initial conditions for this system at the middle of the wall are 

Pl (0) = 0, 7 i(0)=0 



(5i 



The fact that the value of the phase in the middle is not shifted is a trivial corollary 
of the symmetry conditions (|5.3|) . To understand why also p(0) is not shifted, one 
should look at the relation (|5.4j) . Linearizing it with respect to -R(O), p(0), we 
observe that the partial derivative of the left-hand side over -R(O) is zero for the 0— 
order solution while the derivative over p(0) is not. Hence, pi(0) = 0, indeed. When 
the initial conditions ( |5.S| ) are posed, the solutions pi(-z), «i(-z) of the system ( |5.7|) 
should approach zero at z — > 00. We checked numerically that they do. 

Bearing in mind that pi(0) = 0, it is easy to find what i?i(0) is. For z — 0, the 
second equation in (ETBT) is reduced to 



Ri(0) 



9 



m + 



3pH0) 



(5.9) 



where the relation ( |5.5|) and the second equation in Eq. ( |4.3| ) were used. To find the 
shift in the next order, one should expand the equation for /?'(0) further. We obtain 



fttrt fiP^m 3i ?i(°) 4(iV - l)iZo(O) 

m = 6i?2(0) + -rM + — P0T0) — P2(0) 

The second order shift in p(0) can be found from the expansion of Eq. 

3i2o(0)i2?(0) 



(5.10) 



p 2 (o) 



(5.11) 



/3[(0) can be found by differentiating the solution of the linear algebraic system ( |5.6| ) 
at z — 0. Using Eqs. (|4.3| , ^.5|) again and combining everything, we arrive at the final 



(5 
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result 



7](m) = 

(N — l) 3 
162^(0) 



Ro(o) 



1 - 



m + 



3^(0) 



(N-l) 2 
9i? (0) 

m(7N - 1) 



m 



3pr(o) 

32(iV- 1)" 



+ 



3pT(0) 
O[m 3 /Rl(0)} 



For N = 2 when po(0) = p* and -Rq(O) = -R*, the expression simplifies [XT 



rj(m) 



1 - 



2 /4m 



5\ 1/6 



4m 



5\ V3 



9 



0(m 5/2 l 



(5.12) 



(5.13) 



6 Two BPS solutions and Phase Transition in Mass 

Whem m is neither too large nor too small, the Born-Oppenheimer approximation 
does not apply and we are in a position to solve the full system of 4 equations ( |5.1|) 
numerically. We did it for N = 3 and N = 4. As was mentioned in the previous 
section, we use the symmetry relations ( |5.3|) which fix the initial conditions at the 
middle of the wall for the phases and the relation ( |5.4| ) between p(0) and R(0) so 
that only one free parameter is left. We fit it so that the solution would tend to the 
vacuum values ( p.2| ) when z — > oo. 

The situation turned out to be pretty much analogous to what happens at N — 2. 
First of all, the solutions exist only in a limited range of masses. When m is larger 
than some critical value m*, the integral trajectory always misses the vacuum fl5.2|) 
no matter what value for R(0) is chosen. Second, when m < m*, there are not one, 
but two different solutions with a larger and a smaller value of R(0). In Figs. [3], [| 
we plotted the dependence of rj — R(0)/R(oo) on m for both branches. We see that, 
at m = m*, two branches are joined together. This is the reason why no solution 
exists at larger masses. f\ 

The value of the critical mass m* for higher unitary groups turned out to be 
much smaller than that for the SU(2) theory: 



rn 



SU(3) 



.28604 



m 



SU(4) 



.07539 . . . 



(6.1) 



as compared to m* = 4.67059 ... for SU(2). 

Let us discuss now what happens with these two branches in the small mass limit. 
Consider first the upper branch. The illustrative profiles p(z) / p* in the SU(3) theory 
for the upper branch at m — .1 and for m = m* = .28604 (when the solution is 
unique) are plotted in Fig. Comparing the dashed curves in Fig. |5| and Fig. £| we 
see that, for small masses, the solution approaches, as it should, the "Higgs wall" 

4 The presence of two solutions at m < m* and their absence at m > m» can be naturally 



understood by making some simple observations on the phase portrait of the system (5.1) |l 
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ri(m) 




(). 1 



0.2 



Figure 3: The ratio r\ = R(0)/Rq(0) as a function of mass for the SU(3) theory. 
The dashed line describes the analytic result (|5.12 ) valid for small masses. 



ri(m) 




Figure 4: The same as Fig. EJ for the SU(4) theory. 
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solution studied in Sect. 4. The numerical findings for i](m) are in an excellent 
agreement with the approximal analytic result ( |5.12|) both for SU(3) and for S77(4). 

The result Q5.12p allows one to understand why the values of m* fall down with 
N and to make an estimate of m* at large N. Indeed, the critical mass m* and 
the value of the mass where the Born-Oppenheimer expansion in Eq.( p.l2[ ) breaks 



down should be of the same order. For large N, the expansion parameter is k ~ 
-/V 2 (m/-R*) ~ N 2 m 2 / 3 . Assuming k ~ 1, we arrive at the conclusion that m* falls 
down as 

m*(N) oc N~ 3 

in the limit N — ► oo. 

The lower BPS branch is something new which was not and could not be seen 
in the framework of the effective Higgs theory (|2.17f) . The latter was obtained by 
freezing down to zero the logarithm in the full potential ( |2.5| ). But, as was already 
mentioned, we could equally well set = 0. Actually, we see, indeed, that -R(O) — > 
when m — > for the lower branch, i.e. the solution passes in the vicinity of the 
chirally symmetric vacuum on its way from one chirally asymmetric vacuum to 
another. 

A similar lower BPS solution exist also for N = 2, but in that case it is much 
simpler approaching in the limit m->0a combination of two real walls separated at 
a large distance (The phases a(z) and (5{z) are changed, of course, so that the whole 
solution is complex, but the change occurs in the central region where the absolute 
values of the fields R and p are very small.). For N > 3, it cannot be true because 
the energy of a complex BPS wall is not just twice the energy of a complex BPS 
wall, but is lower. One can say that, at small masses, two real BPS walls with the 
energy ( |2.12j) attract each other at large distances and two different bound states 
with the same energy fl2~T3l) are formed. 



7 Lower BPS branch at small masses 

The properties of this new complex BPS branch are rather peculiar. They are 
worth studying in some details. As was also the case for the real walls, this solution 
penetrates a high and narrow potential barrier so that the central region of the wall 
finds itself in the vicinity of the chirally symmetric vacuum. However, in contrast 
to the case N = 2, it always passes the minimum = x — at a finite distance: 
though lim m _>o -R(O) is zero, lim m ^o p(0) is not. This can be immediately seen from 

5 This terminology looks somewhat artificial in 4 dimensions where the walls have infinite energy 
involving the area factor, but one could equaly well discuss a 2 -dimensional Af = 2 supersymmetric 
Wess-Zumino model with the superpotential in which case the walls are interpreted as solitons 
and the wall spectrum is just the spectrum, at least an essential part of it (there may be also some 
breathers). E.g., for N = 3 and to < to*, we have in our disposal 6 real bosonic solitons and 
antisolitons with the mass e r and twice as much complex bosonic solitons and antisolitons with the 
mass V3e r , each soliton possessing a superpartner associated with the fermion zero mode dwelling 
on the wall. One may wonder whether the model is exactly solvable in Zamolodchikov's sense ? 
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Figure 5: The ratio r(z) = p(z)/p* in the SU(3) theory as a function of z for 
171 = 171* (solid line) and for the upper BPS branch at m — .1 (dashed line). 
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the relation ( |5.4j ). For iV > 3, the right-hand side is non-zero, and, when -R(O) — > 0, 
p(0) tends to the value 

/iVcos^ / 4 \ 1 / 2N iNcos^ , . 

p m =o(0) S Po = (^) = ^f* (7.1) 

Let us look first at the numerical solution for p(z) and R(z) for the lower BPS wall 
at iV = 3, m — .01 shown in Figs. [|, 0. We see that the wall consists of three 
distinct regions. First, when going out from the middle point, p(z) goes gradually 
up while R(z) stays practically at zero. Then, in the intermediate region, R(z) 
changes rather abruptly while p(z) is practically not changed. Finally, p(z) and 
R(z) change together in a smooth concerted way until the trajectory levels off at 
the asymmetric vacuum values p* and R*. 

This complicated pattern can be understood and described analytically. The 
analysis is quite parallel to what has been done earlier for the real walls. In the 
limit m — > and while R(z) stays small as is the case in the central wall region, the 
system ( |5.1| ) is greatly simplified and acquires the form 

p> = (JV_ l)mpsin7 . . 

Y = 2(N - l)mcos7 { > 

It has the analytic solution 

p(z) = p cosh 1/2 [2m(Ar-l) 2 ] , . 

sin7(^) = tanh[2m(iV - l)z] 1 ' ' 

Indeed, the numerical solutions for p(z) and 7(2) follow the formulae ( |7.3|) rather 
closely up to the turning point at z ~ 20. On the other hand, for the outer region 
of the wall we are in the Higgs region. The fields R(z) and j3(z) are frozen as is 
dictated by Eq . (|5~5|) , and it is the equation system ( |4.3|) rather than Eq. (|7.2| ) which 
describes the wall dynamics. The change of the regime occurs when the trajectory 
( [7.3|) intercepts the integral trajectory of Eq. ([4.3|) . For > 3, the dependencies 
p(z) and 7(2) in the Higgs phase can be found only numerically, but the trajectory 
in the (p, 7) plane is described analytically according to Eq.( |4.6|) . The interception 
occurs at 



Px = I I . 7x = 7T (7.4) 

/C0S 2(]V3T) 2(iV-l) 

In the transitional region, we should fix instead p(z) and a(z) = 7(2) /2 + 7T / (2N) 
to their values in Eq.([7.4|) after which the system ( |5.1| ) is reduced to 

R' = -2R 2 \n[R 3 pl iN - 1] } (7.5) 
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Figure 6: Profile p(z) of the lower BPS branch at N = 3, m — 0.01 
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Figure 7: Profile R(z) of the lower BPS branch at iV = 3, m — 0.01. 



21 



and (3 is fixed at the value 

/3 X = 



N-l ( 7T 

l 7x + N 



71 

3N 



71 



(7.6) 



After a trivial resclaling, the equation (|7.5|) just coincides with the BPS equation 
( p.2| ) describing the real wall in pure supersymmetric gluodynamics. As was the 
case for the real walls, the transitional region, the region where the potential barrier 



between the Higgs phase 



X 



-2{N- 



x )/ 3 and the symmetric phase = is pene- 



trated, is rather narrow. Its width in z is of order p 2 ^ oc m~( N ~ 1 ^ 3N . On the 

other hand, the characteristic width of the central and of the outer region is large 
oc 1/m. 

Note that, for N = 2 when our wall is a superposition of two distant real walls, 
its structure is simpler. A wide outer region is absent simply because the values 
( [7.4D where the transition occurs coincide in this case with the asymmetric vacuum 
values. But the wide central region and the narrow transitional region are still there. 

The full energy of the configuration described coincides, of course, for small 
masses with the BPS bound ( |2.13j ). One can find out that the fractions of the 
wall energy f c _, /trans, and / ou t carried, correspondingly, by the central region, the 
transitional region, and the outer region are f\ 



fc - 
/out 

For N = 3,f c 



tan 



2(N-1) 



tan 



A? 



1 - fc ~ f 



fti 



TV sin § 



trans. 



•58, /, 



trans. 



.34, and / out « .08. 



(N - 1) cos 



2(N-1) 



N cos jj 



N-l 



(7.7) 



8 Discussion 

Our main result is that, while the real BPS domain walls connecting the chirally 
symmetric and a chirally asymmetric vacua are present at all masses, the complex 
BPS walls interpolating between different asymmetric vacua exist only for small 
enough masses m < m*, m* being given in Eq. ( |6.1| ). A kind of phase transition 
associated with the restructuring of the wall spectrum occurs. [] It makes sense to 
express the result in invariant terms and to trade ^sqcd for an invariant physical 
quantity such as the gluino condensate E = | < Tr A 2 > | in a chirally asymmetric 
vacuum. Restoring the dimensional factor JXgQCD N m Eqs.( [L.2[, |2.9| ) and combining 
it with Eq. (|6.1|) , we obtain 

mf® « .085E 1 / 3 , mfM w 0.033E 1 / 3 (8.1) 

6 To avoid confusion, note that, by derivation, this result is valid only for N > 3. 

7 Needless to say, it is not a phase transition of a habitual thermodynamic variety. In particular, 
the vacuum energy is zero both below and above the phase transition point — supersymmetry is 
never broken here. Hence E vac (m) = is not singular at m = m*. 
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We cannot, however, claim that the quantitative estimates in Eg. flS.ip are, indeed, 
quite correct. The matter is that the particular values of the factors in Eg. (|B.1|) 
are sensitive to the form of the kinetic term in the effective lagrangian which, in 
contrast to the potential term, is not fixed rigidly by symmetry considerations. Just 
multiplying, say, the standard kinetic term for the X superfield C^ in = 1/4/ d 4 9XX 
by a numerical factor would change the particular values ( |8.1| ) of m, ||11|| . The 
effective lagrangian may also involve complicated kinetic structures with higher field 
derivatives oc (<90) 4 etc. 

The effective TVY lagrangian implies the presense of the chirally invariant phase 
0. Recently, this conclusion has been criticized [TJ, |I6| . In particular, it was shown 



in Ref . |L5| that the assumption that the spectrum of the chirally symmetric phase 
in the supersymmetric gluodynamics involves the massless particles associated with 
the colorless composite field $ in the effective lagrangian does not conform with the 
't Hooft anomaly matching conditions of some global discrete anomalies. 

However, this does not exclude the possibility that the chirally symmetric phase 
still exists, but its spectrum has nothing to do with that extracted from the naive 
VY effective lagrangian with the standard kinetic term. Indeed, for the dimension 
of the above mentioned terms with higher derivatives of <fi to be correct, they should 
involve some powers of <ft i n the denominator (no new dimensionful constants can 
be introduced: that would spoil the conformal properties of the lagrangian). The 
presense of such singular terms would modify the spectrum of the chirally symmetric 
phase completely [17]. Actually, the dynamics of the symmetric phase is quite 
unclear by now. This is an extremely interesting problem to solve but, at the 
moment, we do not have insights in this direction. 

We want to emphasize, however, that though higher-derivative singular terms 
may destroy completely a naive spectrum picture, they do not affect the conclusion 
on the existence of the symmetric phase. When the fields are static and homogeneous, 
the form of the effective lagrangian is extracted quite rigidly from the requirement 
that the conformal and chiral anomaly of the original theory are reproduced cor- 
rectly. Also, the dynamics of the domain walls depend largely on the region between 
the vacua where the extra terms with higher derivatives are not singular. Thereby, 
they should not modify the structure of the walls revealed in Refs. [LT], |TT| and in 
this paper. At least, this is our guess and hope. In particular, a qualitative conclu- 
sion on the absence of a smooth transition between the small mass region and the 
large mass region in supersymmetric QCD is correct. 

In this paper, we studied only the BPS solutions, but, as it was done earlier 
for SU(2) |IH , one may and actually should study also non-BPS walls, the field 
configurations which satisfy the equations of motion but not the first-order BPS 
equations. This work is now in progress |]18|| . Our preliminary results display that 
the picture is roughly the same as for SU(2): in some range of masses m* < m < m**, 
a non-BPS complex domain wall presenting a local minimum of the energy functional 
exists. There are also sphaleron wall solutions. For m > m**, the complex walls 
disappear altogether from the spectrum. 
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